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Abstract — In this paper, based on first-principles 
calculations within the pseudopotential theory calculated the 
band states of crystals SnS with rhombic lattice structure. 
According to the calculated band structure of the compound SnS 
valence band can be divided into three subgroups. Analysis of 
wave functions of valence states shows that the lowermost group 
separated from the main group of valence band by a wide energy 
gap about 6eV, arises from s - states of anion. The next group, 
consisting of four bands and located near-7eV is due to s-states 
of Sn. The uppermost group of twelve bands located in the range 
between 0 and - 5eV appears p-states of cation and anion. 
Analysis of origin valence states is suitable with photoelectron 
emission data. The width of forbidden zone computed using 
energy-band structure is suitable to its experimentally 
determined value. Spectral dependence of refraction index, 
reflection coefficient, of the real and imaginary parts of optical 
electric conduction, of the characteristic function of electron 
energy losses and the effective state density of compound for ellc 
and e±c polarizations have been done. 

Index Terms — Energy-band structure, optical properties, 
valence states. 


I. INTRODUCTION 

The present paper contains the results of computing 
energy spectrum and optical functions of SnS compound. The 
computations were made using pseudopotential method 
which is one of the principal for computing energy spectrum 
of charge carriers of semiconductors. This method is 
extensively described in papers [1-6]. 

II. THE EXPERIMENTAL METHOD 

The pseudopotential theory is based on three 
fundamental physical approximations. 

1. The first approximation is the approximation of 
self-matched field. Interaction between electrons in this 
approximation is characterized by a particular average 
potential which itself depends on the states of electrons, the 
mentioned electron states in turn are determined by average 
potential. 

2. In the second approximation all electron states are 
divided into inner shells (“the core”) and of conduction band 
wave functions of the inner shells are supposed to be strongly 
localized. 

3. The third fundamental approximation is the use of 
perturbation theory for electrons in conduction band. 
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Non-local ionic pseudoptentials in configuration space were 
constructed according to a scheme proposed in papers [4, 6]. 

Paper [7] demonstrates symmetry elements of 

D l 2 h (P C mn) spatial group, Brullien zone, tables of 

non-reducible representations of groups of wave vectors of 
symmetrical points and lines of Brullien zone and also 
describes rules of selecting dipole transitions for some 
symmetrical points and lines. 

We have performed the computation of SnS 
compound energy-band structure by pseudopotential method. 
Non - local ionic pseudopotentials in configuration space 
were constructed according to a scheme proposed in paper 
[5]. When computing energy- band structure of this 
compound the compound the screening of ionic charge as well 
as exchange - correlational effects were considered within 
dielectrical formalism in conformity with Hubbard -Sham’s 
model with some selective distribution of charge around each 
ion. About 2500 plane waves were used in expanding wave 
function. In doing so maximum kinetic energy of considered 
energy of plane waves was 20 Ry. 


The parameters of optimized lattice: 


Lattice parameters, 

(A) 

SnS 

Theor. 

Exper 

a 

4.246 

4.334 

b 

10.943 

11.200 

c 

3.922 

3.987 

x (cation) 

0.1155 

0.1198 

y (cation) 

0.1175 

0.1194 

x (anion) 

0.4745 

0.4793 

y (anion) 

0.8535 

0.8508 


As seen from Fig.l, valence band of SnS consists of three 
groups of bands. The characteristic pair wise arrangement of 
valence bands is consequence of crystal lamination. Analysis 
of wave functions of valence states shows that the lowermost 
group separated from the main group of valence band by a 
wide energy gap about 6eV, arises from y - states of anion. 
The next group, consisting of four bands and located near - 
7eV is due to y-states of Sn. The uppermost group of twelve 
bands located in the range between 0 and - 5eV appears 
because of p-states of cation and anion. Analysis of origin of 
valence states agrees well with photoelectron emission data. 
The width of forbidden zone computed using energy-band 
structure is in good agreement with its experimentally 
determined value. 

It is known that optical properties of SnS are not 
sufficiently investigated experimentally. For this reason 
theoretical investigation of several optical properties of the 
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said compound is needed. With this object in view we real 
parts have conducted computations of imaginary and real 
parts of complex dielectric penetrability in the course of 
polarizations of parallel and perpendicular optical axis and 
have computed other optical functions as 
well. 

SnS atom 0 size 0.20 
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Fig, 1, Energy- band structure of 

compound 


where CD n - plasma frequency for electrons; n- average 

density of electrons in crystal. The computation shows that 
fl(D p equals to 14,25 eV. 

The histogram was plotted with a step of ^0,2eV. 
All transitions v^c with energy up to 15eV had been 
considered. Starting from small structure in the histogram 

about 15eV the £, (&>) dependence was extrapolated by the 

known formula: 

^ i »oo ~ 3 • 

CD 

The real part of dielectric penetrability is computed 
from integral dispersive Kramers-Kroning relation 

/ \ , 2 r . ( \ dco 

£ r {co) = \ + — P\co sAcq )•— - 

71 J Q CD - CO 

symbol P here stands for integral in a sense of the main value. 

The effective number of valence electrons per atom 
participating in transitions with the energy of E <Fl CD is 
determined as: 


n 


eff 


m e 

O 2 2 

2n e 




N 


cell 


V. 


cell 


9 


where n a -density of atoms in crystal; e - elementary charge; 


III. RESULTS AND DISCUSSIONS 


When determining frequency dependence of the imaginary 
part of complex dielectric penetrability we have employed 
relation [8]: 


a 2 2 


m CD 


vc zB 


2 d 3 k 

{2xf 


M CT (k) 


S(E c (k)-E v (k)-tuo), 



here e • 


M v ( k ) = e - J y/* ck (r X- ihV }// vk (r )d 3 r, where 


an integral in the right part denotes matrix element of pulse 
operator p = —itlSZ ; indices v and c numerate states of 


valence band and conduction one respectively; k- wave 
vector; e - unit vector of polarization. Integration is made 
with respect to the volume of crystal elementary cell. 

In (1) we have replaced integration with respect to BZ by 
summation in elementary cell of inverse lattice. In this process 
the elementary cell was divided into 8 part of equal volume 
and k points were selected in these parts randomly. 1280 
points were taken in all as a result of which 
a smooth histogram was obtained. Then this histogram was 
normalized by the following formula 


I cos i (cojdco = 

0 


n 2 
co n 

2 ‘ 


n 47rn e e 2 
2 m e 


m e - electron mass. 



Fig. 2. Spectral dependence of the imaginary (ci) and real (c r ) 
parts of dielectrical penetrability of SnS compound for e||c 

(a) and e_Lc (b) polarizations 
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The results of computing the above mentioned 
optical constants in the 0-8 eV energy range are displayed in 
Figs. 2-5. As follows from Fig. 2a, a maximum of the main 
peak in s r (co) spectrum at e||c polarization is observed when 
the energy is 2,20 eV, when the energy is 3,37 eV a weak 
maximum is noted. c r (E) goes down to the energy of 6,46 eV 
after which increase takes place. 



-1 0 1 2 3 \ 5 E J 3 



Fig . 3. En ergy dependence of refra ctio n in dex (a) and 
reflection coefficient (b) of SfiS corn pound at e/Jt and 


eJ.c polarizations 

At 4,25 eV c r (E) inversion occurs. At e_Lc polarization c r (E) 
increases to the energy of 2,23 eV , reaching its maximum 
value (24,4) drops to 3,13 eV, a small growth and decrease in 
c r (E) are observed. Fig. 2b demonstrates dependencies of the 
imaginary part of dielectric penetrability c r (E) on energy at 
e||c and e_Lc polarizations. 

As seen from Fig. 2b at e_Lc polarization Ci(E) at first 
rises to 13 when the energy is 3,4 eV, attaining a maximum 
and falls to 7,97 at 6,14 eV . Further, at the energy of 6,23 eV 
two maxima in c r (E) dependencies at e||c polarization happen 
at 3,1 and 4,6 eV and a minimum is observed at 3,54eV with a 
decrease in Ci(E)down to 6,46. 

As follows from Fig. 2b, at e||c polarization two 
clearly marked maxima are seen on Ci(E) curve at 3, 1 and 4,66 


eV and a minimum occurs at 3,54 eV. After 4.66 eV a 
decrease in 8i takes place, Si reaches its 
maximum value of 13,6 when the energy is 3,4 eV. At e_Lc 
polarization two maxima are evident on Ci(E) curve at 3,4 and 
6,8 eV and deep minimum at the energy of 6,14 eV is 
witnessed. On the whole, at both polarizations changes 
according to energy are identical. 

The imaginary part of refraction index is determined 
by the following formula: 

k = J-(-£ r + 4sl + £•) 

2 

while reflection index at normal falling of light is determined 
by the formula below: 

R = ( ft - 1 ) 2 +£ 2 

(n + lf+k 2 ' 

Fig. 3 gives the results of computing refraction (k) and 
reflection ® indices at e||c and e_Lc polarizations. As follows 
from Fig. 3a, at e||c polarization k rises from zero to 2,8 and at 
the energy of 3,32 eV a slightly marked maximum and a weak 
decrease are observed. 




Fig. 4. Spectral dependencies of the real c^fE} {a} 
and imaginary parts cr{E) (bi of optical electric 
conduction of SnS compound for e/Jb and eJ.c 
polarizations 
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At e±c polarization two maxima are present on k(E) 
curve at the energies of 4,79 and 6,55 eV and two minima at 
the energies of 3,51 and 5.97 eV, at e_Lc polarization an 
increase in k occurs from zero to 3,5 with differentials. 
Reflection index at e||c polarization has two weakly marked 
maxima at 3,25 and 5,1 eV (Fig3b) 

On the whole, R change with energy is little-from 0,35 
to 0,5. At e_Lc polarization three maxima are witnessed on R(E) 
curve at 2,17; 4,79 and 6,52 eV and two minima at 3,46 and 
5,9eV. And in this case R change with the growth of energy 
takes place I the 3, 5 -5, 5 interval. 

The real and imaginary parts of optical electric 
conduction of SnS are found with the help of the following 
formulae: 



CO£ i _ COS r 

\n 4 n 


The results of computing the real and imaginary parts of 
optical electric conduction are given in Fig.4. As seen from 
Fig.4b, two weakly marked maxima are present on a r (E) 
curve at e_Lc polarization when the energies are 3,28 and 6,69 
eV and minimum at 5,9eV a r goes up in the 0-3,28 eV energy 
range, it remains practically constant in the 3,28-5,9 eV range, 
its drastic rise from 3,5 to 8,5 is observed in a narrow range 
from 5,9 to 6,7 eV, further a weak decrease occurs. 

At e||c polarization a structure is found on a r (E) curve at 
3,02 eV, a sharply marked maximum is observed at 4,59 eV 
and a weak maximum is present 6,57 eV. The a r (E) changes at 
e||c polarization in the 0-8 eV energy interval take place from 
zero to 6. Fig.4a shows Gi(E) dependencies at e||c and e_Lc 
polarizations. As follows from Fig.4a, at e_Lc polarization two 
sharply pronounced maxima are observed on the CTi(E) curve 
at 2,23 and 3,5 eV and two minima at 3,01 and 6,41 eV. On 
the whole, structures and electric conduction intervals are 
present on the ai(E) curve. A change in in the 0-8eV 
interval ranges from -2,5 to 4,5. At e||c polarization two 
maxima are registered on the Gi(E) curve at 2,23 and 3,56 eV, 
a minimum is revealed at the energy of 3,3 eV. A change in Oi 
at e 1 1 c polarization in the energy interval under investigation 
ranges between -2 and 4. 



cos i _ cos r 
4n 4 n 


The characteristic function of electron energy losses is 
determined in this way: 



(-) 

UJ 



sf +£ 


2 • 

r 


The spectral characteristic of the imaginary part of 
the inverse value of complex dielectric penetrability 

(/.(-«-■)) is displayed in Fig. 5 a. It is seen from here that 

no clearly marked maxima and minima are observed at e||c 

polarization, i.e (/.(-«-■)) rises over a wide range of 

energies from 0 to 6,3 eV and decreases over a narrow energy 
range. 

At e_Lc polarization {i m (— S 1 )) at first grows, an 
arm is detected in the 3,0 - 4,7 eV interval, a clearly marked 


maximum is present at 6,01 eV a minimum at the energy of 
6,75 eV. 


The effective state density is expressed by the 
formula below 

8eff = AM 2 ’ 

and the results of computing at e||c and e_Lc polarizations are 
given in Fig. 5b. As follows from Fig. 5b, at e||c polarization 

g e fj slightly increases in the 5.82-6.83 eV range, smeared 

maxima are witnessed at the energies of 3,82 and 5,36 eV and 
minimum at 5,82 eV. An arm is observed on the g e ^ versus 

- energy dependence curve in the 3-3,4 eV range at e_Lc 
polarization, a maximum takes place at 4,67 eV and a 
minimum occurs at 5.95eV. An abrupt rise in the 6 -6,5eV 

range and a slight decrease in g e ^ are noted. 


IV. Conclusion 

So, in this report the spectra of optical fundamental functions 
of the semiconductor SnS have been calculated for two 
polarizations ellc and e_Lc in 0 - 8 eV energy range for the first 
time. Parameters have been defined, effective masses of 
electrons and holes of compound were calculated on the 
basis of the band structure. The results obtained allow a 
detailed analysis of the optical properties and electronic 
structure of SnS 
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